Decay process of quantum open system at finite-temperature 
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Starting from the formal solution to the Heisenberg equation, we revisit an universal model for 
a quantum open system with a harmonic oscillator linearly coupled to a boson bath. The analysis 
of the decay process for a Fock state and a coherent state demonstrate that this method is very 
useful in dealing with the problems in decay process of the open system. For finite temperature, the 
calculations of the reduced density matrix and the mean excitation number for the open system show 
that an initial coherent state will evolve into a temperature-dependant coherent state after tracing 
over the bath variables. Also in short-time limit, a temperature-dependant effective Hamiltonian 
for the open system characterizes the decay process of the open system. 

PACS numbers: 05.40.-a, 03.65.Yz, 42.50.-p 



I. INTRODUCTION 

Actually, no quantum system can be isolated from 
its surrounding environment completely Q . When it 
is coupled to the environment, the system (also named 
open system) energy dissipates into the environment ir- 
reversibly. There exists two approaches to deal with the 
decay process of the open system, i.e., the system plus 
bath model 3] and effective Hamiltonian model In 
a simple model of a harmonic oscillator interacting with a 
bath of many harmonic oscillators, the relation between 
these two approaches is clarified by Yu and Sun [1, @] ■ 

In quantum optics and macroscopic quantum phenom- 
ena, the mesoscopic or macroscopic states can be rep- 
resented by coherent states. Due to the quantum fluc- 
tuation and thermal fluctuation of the bath, it is very 
important to study how the open system evolves with 
an initial Fock state Q and an initial coherent state 
especially at finite temperature. 

Usually, the master equation method is adopted to 
treat the decay process of the open system 10], time evo- 
lution of the reduced matrix for the open system can be 
solved analytically and numerically. In this paper, start- 
ing from the formal solution to the Heisenberg equation, 
we can construct the reduced density matrix for the open 
system, this method has been applied to deal with some 
problems in quantum open system [(| 0, HTI - flij ]. Here, 
an exactly solvable model of one boson interacting with 
a bath of many bosons in rotating-wave approximation 
is used to analyze the decay process of the open system. 

The paper is organized as follows. In sec. II, we describe 
an universal model for quantum open system, a single 
mode boson interacting with a bath of many bosons 
through system-bath coupling. In sec. Ill, we revisit the 
decay process of a Fock state and a coherent state at zero 
temperature. In sec. IV, we study the finite-temperature 
decay process of the open system. Finally, we summarize 
our results. 



II. SYSTEM-BATH MODEL 

The quantum decay process is characterized with the 
vanishing of diagonal elements of the time evolution of 
reduced density matrix for an open system. The simplest 
model is adopted that the composite system consists of 
the system (a single mode boson) interacting with the 
bath of multi-mode bosons. Using the method in Refs. 0, 
0, ITll - flij , we can construct the reduced density matrix 
and calculate the mean excitation number for the open 
system. 

The model Hamiltonian for this composite system is 



H = hw h b% - 



where b (a,j) and b^ (aj) are bosonic annihilation and 
creation operators for the system (bath), and the £j are 
complex coupling constants parameterized via the spec- 
tral density function J(w), i.e., 



j m =y,ioi 5(ui-wj). 



(2) 



Here we have applied rotating-wave approximation on 
the model Hamiltonian H . 

Now we can use the time evolution of operator in 
Heisenberg picture to construct a wavefunction (or a re- 
duced density matrix). These well-known solutions of 
Heisenberg operators are given in Ref. [ll[ 

b{t) = t* (t) 6 (0) + 5>i(t) %(()), 
3 

a,j(t) = e'^aj (0) + Uj (t) b (0) + ^ Vj, s (t) a 8 (0) . 



Where some coefficients in the above equation are 



u(t) 



v, (t) = fce"*" 



e 2 e 



LOl, — UJj — l- 
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Uj (t) = c; e -^ 4 



(3) 
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When the spectral function of the bath is assumed as 

■/(«) = J, (4) 
the calculations show that 

X)l«j(*)| 2 = l-e-^. (5) 



and |u (i)| 2 + \ vj (t)\ 2 = 1 which maintains the com- 
mutation relation [b (t) , (£)] = 1. 



Here we have taked into account the relation 



(9) 



where A, = u*(i)aj(0). According to the result 
in Eq.®, we obtain the probability in the state \n), 
P n (t) — e - ™ 7 ' (seen in Ref. ||). Here the decay time 
of the open system is defines by 



T = 717- 



(10) 



III. ZERO-TEMPERATURE DECAY PROCESS 



In quantum optics, the study of Fock state and coher- 
ent state is very important. At zero temperature, the 
decay of the Fock state and coherent state is presented 
with master equation method in Refs. @, Q. Here, using 
the method in Refs. @, 0, we revisit the decay of 

a Fock state and a coherent state. 



A. Decay of a Fock state 



The initial state of the composite system is assumed 



|V(0)) = |n)®|{0 j }>, 



(6) 



where n stands for the excitation number for the Fock 
state \n) and the bath is in vacuum state |{0j}). The 
ground state of the composite system is |0)<g)|{0,,}) , which 
is invariant under the operation of Unitary evolution op- 
erator U(t), i.e., U(t) |0)®|{0j}) = |0)®|{0j}). Then we 
can construct the time evolution of the composite system, 



IV (*)> = U(t)\i>(0)) 



|0) ® |{0,}) 



Where the operator B (t) is defined as B (t) — 
U (t) b (0) W (t). Using the method in Ref. [3, we can 
replace i with —i in Eq.Q and get 



fit (t)=„(*)6t (0) + X;«j(*)oJ(0). 



(7) 



Tracing over the bath variables in density matrix for 
the composite system p(t) = \ip(i)) (ip(t)\, the reduced 
density matrix for the open system is calculated as 

71 

Ps (t) = p m (*) \m) (m| 
m=0 

where the probability in the state \m) for m — 0...n is 
n! 



Pm{t) 



(n — m)!m! 



-m^yt 



1 - e 



(8) 



B. Decay of a coherent state 

Similar as the decay of a Fock state, for an initial co- 
herent state | a), we construct the time evolution of wave- 
function for the composite system, 



hK*)> 



_ p aB^t)-a*B(t) 



|0) ® |{0,}) 



.(*)>®ni a ^(*)) 



Then the corresponding reduced density matrix is 

p s (t) = \au (t)) (ou (t)|. 

Through the calculation of the mean excitation number 
for the open system, 

N = Tr (p. (t) 6+6) - \a\ 2 e-~<\ 

the decay time of a coherent state \a) is obtained (seen 
in Ref.[|), 



T 1 = 7. 



(11) 



When the bath is in excited state (coherent state 
|{Aj})), the master equation method is not valid. Us- 
ing our method, we can easily obtain the time evolution 
of wavefunction for the composite, 

|V (*)> = e aB *- a * B \0) ®~[[e^ A >- x > A i \0j) 

3 

= ® |{Mj (*)}> • (12) 

where the two coefficients in above equation are 

p (t) = au (t) + XjUj (t) , 
i 
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IV. FINITE-TEMPERATURE DECAY PROCESS 

In previous section, it shows that when the bath is 
prepared in a vacuum state, the quantum vacuum fluc- 
tuation of the bath will induce decay process of the open 
system. Here we will demonstrate the finite-temperature 
influence on decay process of the open system when the 
bath is in a thermal equilibrium state. 

Initially, the density operator for the composite system 
is written as a direct product 

p(o) = (|V(o)> (V (o)|)®/»b. 

Here the initial state of the open system is assumed in a 
coherent state, \ip (0)) = \a). In coherent-state represen- 
tation, the bath at thermal equilibrium is described by 
the density operator 



PB 



where the mean excitation number in the j'-th mode of 
the bath with the frequency uij is 



= (e^ - 1) 



(13) 



Using the result in Eq.(|12j). the reduced density matrix 
for the open system is calculated, 



p s (*)=n 



d 2 X 



Im(*)> <m(*)I 



(14) 



It can be used to characterize the decay process of the 
open system. Then the probability of the open system in 
the initial state \tp (0)} = |a), i.e., diagonal elements of 
the reduced density matrix, can be obtained 

(^(G)|p s (i)hK0)> = <a|(|V. (t)) (t)\)\a). 

Here the temperature-dependant wavefunction for the 
open system is obtained, 

(*)> = $ (T, tyi \a ((« (t) - 1) $ (T, t)~i + l) ) . 

(15) 

It shows that at finite temperature the state of the sys- 
tem evolves into a coherent state \ip s (t)) according to 
a coherent state (\a)) initially. In Eq. (fT5|) . the time- 
(temperature-)dependant term is denoted as 



$(T,t) = l + J2nj\ Uj (ttf 



(16) 



Changing the sum into an integral in the above equation, 
we have 



2tt 



dio 



R -3£—i(u b -w)t _ i 
K-c) 2 + (i) 2 



Here the term (e^ u — l) 1 is a slow-varying function of 
the frequency ui, thus we can take it out of the integral 
and get 

$(T,t) = l+nth(l-e-T*) , 

where the mean excitation number with the frequency ujb 

is (c" ftu » -l)~\ 

To demonstrate the decay process of the open system, 
we now calculate the mean excitation number for the 
open system. A decay process means that the mean ex- 
citation number reduces with the time increasing. Apply- 
ing the time evolution of the wavefunction for the open 
system in Ea. dT5|) . we can calculate the mean excitation 
number to characterize the decay process of the open 
system, i.e., 

N = (V.(t)|6 + 6|V-(*)) 

For low temperature (/? — > oo), we approximately have 
W[h — e~P hi ^ b . Then the mean excitation number of the 
open system becomes N ~ |a| 2 e~ 7 *. It shows that, in the 
low-temperature limit, the above mean excitation num- 
ber N exponentially decays as the time t increases and 
approximately does not depend the temperature T. That 
is to say, the vacuum fluctuations dominate the decay 
process of the open system. 

For high temperature (/? — » 0), we have nth — 
(fihcob) 1 and then the excitation number becomes 

A~/3^ fc |a| 2 (l-e-^) _1 oc I. 

It shows that, in the high-temperature limit, the rise of 
the temperature will accelerate the decay of the open 
system. 

Assume the system initially in a state \ip s (0)), it can 
be expanded into the superposition of a series of coherent 
states, i.e., 



Ms (0)) = / — C a \a) , 

7T 



(17) 



where the amplitude of the probability in the state \a) is 
C a = (a\ip a (0)). Then the corresponding reduced den- 
sity matrix will be 

Ps (*) = / |^. Q (*)> J ^c; (t)| , (is) 

where for £ = a, /?, we have defined 

(*)) = $ (T, |$ ((« (t) - 1) $ (T, t)~* + l)) . 

In short-time limit, jt -C 1 and $(T,i) ~ e^ 7 *, the 
temperature-dependant wavefunction of the open system 
becomes 



(*)) 



-iH B jft 



10- 



(19) 
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According to the results in Eq. (|18ll9l) . the corresponding 
temperature-dependant reduced density matrix for the 
open system in Eq.([T8|) can be rewritten into 



Ps (t) 



(0)) <V» (0)1 ■ 



where the temperature-dependant effective Hamiltonian 
for the open system is obtained, 



H, 



eff 



- i— J b^b — i-nthl- 



(20) 



Consider the the open system initially in a Fock state, 
IV's (0)) = \n), the corresponding time evolution will be 



(*)) 



and the mean excitation number becomes 

N ~ ne- (Kir+ ™ )7 * (21) 
in which the decay time of the open system is 

t- 1 = {WK+n) 1 . (22) 

In a same fashion, consider the open system initially 
in a coherent state, \tp s (0)) = \a), the open system will 
evolve into a temperature-dependant coherent state 



|Va (*)) - e-* 11 ^ \a) 



oe 



Then the corresponding mean excitation number is cal- 
culated as 



N = |a| e 
with the decay time 

= (Jhh + 1)7- 



(23) 



(24) 



The results in Ea.([22|) and Eq.(|24)l show that thermal 
mean excitation number of the bath nth accelerates the 
decay of the open system. When at zero temperature, 
ritjl = 0, the results in Eg. (f2"2")l and Eg . (|2"4")) come back to 
the results at zero temperature in Eo. dTUl) and Eg.tlTTI). 



CONCLUSIONS 



In this paper, we have studied an universal model for a 
quantum open system with system- bath coupling. Using 
the method in Refs. [1, 0, EH3 ' not master equation 
method, we revisit the decay process of the open sys- 
tem at zero temperature. Specially in short-time limit, 
a temperature-dependant effective Hamiltonian for the 
open system is obtained (seen in Eq. (|20"]) ). For finite 
temperature, it can characterize the decay process of 
the open system very well. The decay time shows that 
the temperature T will accelerate the decay process of 
the open system. In addition, a temperature-dependant 
wavefunction for the open system is easily obtained (seen 
in Eg . (TT51) ) . it shows that at finite temperature, the open 
system will evolve into a temperature-dependant coher- 
ent state from an initial coherent state. In summary, this 
method in Refs. @, 0, ITiTflij is very useful to deal with 
the problems of quantum open system. 
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